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The DO Collaboration has reported an excess of roughly one percent of \x~ \x~ 
pairs over /i + ji + pairs in pp collisions at a center-of-mass energy y/s = 1.96 

ps| | GeV at the Fermilab Tevatron, when known backgrounds are subtracted. 

This excess, if ascribed to CP violation in meson-antimeson mixing of non- 

& ■ strange or strange neutral B mesons, is about 40 times that expected in the 

Standard Model (SM). We propose a null test, based on a tight restriction 
on the muon impact parameter b, to confirm that this excess is indeed due 
to B mesons. If the asymmetry is due to anomalous CP violation in B s -B s 
mixing then a tight restriction on b would increase by a factor two the net 
asymmetry from neutral B mixing, while the sample of dimuons from neutral 
B decays will be reduced significantly relative to background events. 

<n: 

O ■ PACS codes: 12.15.Hh, 12.15. Ji, 13.25.Hw, 14.40.Nd 

<*; 

I Introduction 

The present understanding of CP violation, based on phases in the Cabibbo-Kobayashi- 
Maskawa (CKM) matrix, predicts a very small asymmetry in the yield of same-sign muon 
pairs due to bb production followed by oscillation of a neutral non-strange or strange B 
meson into its antiparticle. Defining 

sl ~ N++ + N— ' 
the DO Collaboration [IJ[2] has reported an asymmetry of 

A b sl = -0.00957 ± 0.00251 (stat) ± 0.00146 (syst) , (2) 

about a factor of 40 larger than the SM prediction [3] of 

A b sl = (-2.3lg.-i) x 10- 4 . (3) 

This result is interpreted as evidence at the 3.2<r level for anomalous CP violation in the 
mixing of neutral B mesons. 
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The DO analysis employs a large number of systematic checks to guarantee the sta- 
bility of their result. These 16 checks lead to values of A b sl consistent with the nominal 
result (J2J). All these checks used subsets of data involving about equal fractions of back- 
ground events (-Fbkg) from non-genuine bb pair production. To extend such tests, we 
suggest a measurement, based in part on a check already performed by DO, to determine 
if a substantially smaller asymmetry is obtained in a sample depleted in bb pairs. The 
method employs a reduction of the maximum allowed impact parameter of muon tracks 
which should reduce the signal more than it reduces sources of background. Basically, 
we are asking if the DO experiment can infer zero signal from b decays when no signal is 
expected. 

In Sec. II we review kinematics of muons in B decays and show the effect of changing 
impact parameter selection criteria. The relation between impact parameter and DO's 
selection criteria is studied in Sec. III. Implications of strict impact parameter cuts are 
discussed in Sec. IV while Sec. V concludes. 



II Production and decay kinematics 

In pp collisions at the Fermilab Tevatron, B mesons are produced with a distribution 
in energy E such that /?7 = 0(1), where 7 = l/i/l — (3 2 = E/m,B (see, e.g., Figs. 2.2 
and 2.4 of Ref. jl]). The pseudorapidity r\ = — logtan(#p/2), where 9p is the angle with 
respect to the proton beam, is typically less than 0(3) in magnitude. 

The average path length traversed by a decaying B meson in its rest frame is about 
ctb = 0.45 mm, where we take B° and B s lifetimes to be r^ = 1.5 ps and neglect 
differences between the two. Let t denote the B proper decay time in any event. Then 
the normalized distribution in t is w(t) = (1/r) exp(— t/r), satisfying J °° w(t)dt = 1. 
The proper path length ct becomes £ = ft'jct in the laboratory frame. 

Define an axis z* in the B rest frame pointing in the direction of the boost (with 
velocity f3) from the laboratory frame to the B rest frame. Let a muon in B semileptonic 
decay be emitted with an angle 9* with respect to this axis. Since the B is spinless, the 
distribution of the muon will be isotropic in cos 9*. We seek the transformation from this 
angle to the angle 9 in the laboratory frame between the B and muon directions. The 
impact parameter of the muon track with respect to the primary vertex is then ^sin^, 
and its distribution is easily calculated for any value of (3. 

Define the x* axis in the B rest frame to be coplanar with the B boost direction (the 
z* axis) and the muon direction. Then the Lorentz transformation between momenta in 
the laboratory frame (unstarred) and the B rest frame (starred) may be written 

p* = pI ; Pz = i(p: + (3e*) ; E = j{E* + Pp* z ) . (4) 

Here E* = (p x 2 + p* 2 ) 1 ^ 2 , E = (p 2 x + p 2 ) 1 ^ 2 ] we have neglected the muon mass. Then 
since p* x = E* sin 9*, p* z = E* cos 6*, p x = E sin 6, p z = E cos 9, we have 

sin 9* 

srn $ = ~~ r-, 7, ttt • (5) 

7(l + /3cos#*) v ' 

This reduces, as it should to sin 9 = sm9* for f3 = 0. Note that this relation does not 



CD 

. i — i 

m 



1.0 



0.8 



0.6 



0.4 



0.2 



0.0 



1 




i> 


"Is 




-ssl l l l l l l l l 


h 


/ / 




s. 


N 


^v ' ' 


.1 


/ / 










If / 










\ N^ 


1 / 












I / / 












T ' / 










\ \ _ 


W 

' 1 / 










\ \ _ 


-!i/ 










\ \ _ 


i 
i'J 










\ \ _ 
\ \ 

\ \ 


H 










\ \ - 


7 ' 
1 ' 

4_\ 










\ \ - 
\ \ 

\ \ - 
\ \ 

\ \ - 
\ \ 

\ \ 


1 ' 












■ 








\ \ ~ 




. 








\ \ _ 




\ 








\ \ 


■ 


v . 








\ \ - 

\ \ 


- 




'x 






\ \ ~ 


— 






"* . 




\ \ 


; 












- 










"" ' ~ ■ - - . _ \T 


i 


i 


1 


1 




1 1 1 1 1 1 1 1 1 1 1 1 f ' -s 



1.0 



-0.5 



0.0 0.5 



1.0 



cos 



Figure 1: Transformation between cos#* and sin# for 0^y = 0.7 [(3 = 0.573] (solid), 
/3 7 = 1 \J3 = 0.707] (dashed), and fa = 10 [f3 = 0.995] (dot-dashed). Note that 6 -> vr/2 
when cos#* = — j3. 



depend on the muon energy (in the limit that m M may be neglected). We illustrate the 
transformation for several illustrative values of 7/3 in Fig. [TJ 

Using this transformation and the isotropy of muon emission in the variable cos#*, 
one may calculate the average value of sin 9 and hence the average impact parameter as 
functions of ^7: 

(b) =-?(3 (sin 6)cr , 



where 



(sin 1 



sin 2 6* d9* 



7T 



1 



7(l + /3cos#*) 2 1 + 7 



(6) 



(7) 



The result is shown in Fig. [2j 

Specifically, for 7/? = (0.5, 1,2) one has (b) = (167,293,437) /xm. Fitting by eye the 
exponential tail of an impact parameter distribution studied by the CDF Collaboration 
(see Fig. 6 of Ref. [5]) we estimate (b) = 350 /zm, which lies in this range. The impact 
parameter for a given 7/? will be distributed as W(b) = (1/(6)) exp(— 6/(6)), where 
/ °° W{b)db = 1. If one excludes values 6 > 6 , the fraction excluded for muons from 
B decays will be J^ W{b)db = exp(— 6 /(6)). Thus, the remaining fraction of dimuons 
from the two B meson decays will be [1 — exp(— 6 /(6))] 2 . We show in Table fl] the 
remaining fraction of dimuon events for various values of (6) and 60 ■ Thus, demanding 




Figure 2: Average impact parameter (b) as a function of 7/3 for a B decay with assumed 
ct = 450 /im. 

Table I: Remaining fraction of dimuon events for various values of average impact pa- 
rameter (b) when events with impact parameter exceeding a value of b are discarded. 
Units are /im. 
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an impact parameter of less than 100 /im should be enough to significantly reduce the 
dimuon signal of two B decays in Refs. [TJ[2] • In the next section we discuss the relation 
of this criterion to those employed in the DO analysis. 



Ill Relation to criteria employed by DO 

The DO Collaboration quotes a precision for primary vertex reconstruction of 20 /im 
in the transverse plane and 40 /zm along the beam direction. Requiring an impact 
parameter of less than 100 /im, as proposed above, thus should not reduce the relative 
fraction of signals from tracks originating from the primary vertex, including kaons, 
pions, and (anti)protons misidentified as muons [lj. However, the employed criteria 
distinguish between the transverse impact parameter relative to the closest primary 



vertex (which we shall call b±) and the longitudinal distance from the point of closest 
approach to this vertex (which we shall call b\\). The DO selection criteria for these 
variables are [TJ[2] 

b± < 3000 /jjD. , 6|| < 5000 fjm . (8) 

We seek the relation between these quantities and the impact parameter b. 

Laboratory coordinates (x, y, z) are defined such that the detector axis (very close 
to the beam axis) is z, the plane of the accelerator is x, and the orthogonal (vertical) 
direction is y [6]. The vertex position (x ,y ,z ) is close, but not necessarily equal, to 
x = y = z = 0. Define another set of coordinates (x', y', z') for each track such that the 
transverse point of closest approach (PCA) has coordinates x' = b±,y' = z' = 0, while 
the vertex is defined to have x' = y' = 0, z' = b\\. 

By definition of the PCA, the linearized approximation to the muon trajectory near 
the PCA lies in the y'—z' plane, making an angle ip with the z' axis. The transverse and 
longitudinal components of muon momentum in the laboratory frame, p± and p z , may 
be used to determine ip: 

P± = V 1 * 1 si n ip i Pif — P 11 cos "^ • (9) 

Let the distance along the muon track from the PCA be denoted by s. The muon 
track's coordinates in the primed system are 

y' — ssinip , z' = scosip . (10) 
The distance d of a point along the muon trajectory from the vertex satisfies the relation 

d 2 = b 2 ± + (s sin tp) 2 + (s cos ^ - b {l ) 2 . (11) 
We seek the minimum value of this quantity with respect to s: 

= d(d 2 ) /ds = 2s -26„ cos ip , (12) 
implying s = b\\ cos ip and hence 

^ min = [6i + (6|| S in^) 2 ] 1/2 = & • (13) 



When ip = 0, b = b±, while when ijj = n/2, b = Jb"]_ + 6|, as expected. It is then a 
simple matter to construct the impact parameter b using the DO track-selection criteria. 

IV Implication of impact parameter cuts 

Indeed one of the 16 analysis variations mentioned above ("Test D" [TJ) changes the 
values of the maximum transverse and longitudinal impact parameters to values tighter 
than (JHJ: 

b ± < 500 /jm , b\\ < 500 /mi . (14) 

However, since one always has b > bj_, one can see from Fig. [2] and Table [I] that this 
is not expected to reduce the signal very much, and it doesn't [TJ. While the tighter 



cuts reduce the number of like-sign dimuon events by a factor of two, the fraction of 
estimated background events is unchanged within errors. (See Table XIV and Eq. (37) 
in Ref. pQ.) 

The much tighter restriction b < 100 /zm should have a noticeable effect, reducing 
significantly the dimuon signal from two muonic B decays relative to background events. 
The remaining fraction of dimuon signal was calculated in Table I to be 8 and 4 percent 
for (b) = 300 and 450 microns, respectively. As a second-order effect, we note that a 
tight restriction on b implies that the charge asymmetry from the remaining neutral B 
decays will be dominated by B s mesons which involve a much higher oscillation frequency 
than B°. If the anomalous charge asymmetry observed in Ref. [TJ[2] originates in B s -B s 
mixing rather than in B°-B° mixing then the negative asymmetry from the net signal 
events would increase by about a factor of two. 

In order to demonstrate this effect, let us take for instance r\> = bo/(b) = 1/3. The 
net asymmetry from neutral B meson mixing is given by 

A b _ fdZ d Ad f s Z s 

Jd^d + Js^s Jd^d + Js^s 

where A d sl and A s sl are asymmetries in B°-B° and B s -B s mixing, respectively. We will 
take B° and B s production fractions fd = 0.323, f s = 0.118 [7J, noting that our result 
will not depend on precise values of these parameters. Neglecting small width differences 
(Ar ds /2r) 2 <C 1, the parameters Z q (q = d, s) are given by 



I X 1111/ /"; i ■ ( '( IS, I ■/" i ' : i 

(16) 



C x q sm(x q rb) — cos(x q rb 



7 = 2 e -n 

q l + xl 



1 + 



l + x 2 q 



where x q = Am q /T. Ignoring small errors in neutral B mass differences, Xd = 0.774, x s = 
26.2 [7], one finds Z d = 0.0029, Z s = 0.264 implying A b al = 0.03^ + 0.97^. This result 
should be compared with A b sl = 0.5lA al + 0A9A s sl [1\ (also ignoring errors) obtained with 
loose restrictions on the muon impact parameter (i. e., r& 3> 1). Thus, a tight restriction 
on the impact parameter would increase by a factor two the net asymmetry from neutral 
B mixing if the asymmetry originates in the B s system. 

A tight restriction on impact parameter eliminates also a large fraction of non-genuine 
semileptonic B decays. This includes hadrons in bottom and charm decays, where the 
hadrons (kaons, pions and protons) are misidentified as muons, and muons from se- 
quential decays of kaons and pions. In order to provide a null test, or alternatively to 
demonstrate a nonzero asymmetry from sources other than B°-B° and B s -B s mixing, 
the number of remaining events, originating mainly in promptly produced hadrons fak- 
ing muons, must be sufficiently large. We note that a contribution to the same-sign 
dimuon asymmetry from kaons, dominating over those from pions and protons, must 
be subtracted from the raw asymmetry in order to obtain the combined CP asymme- 
try in B°-B° and B s -B s mixing. A positive asymmetry of five percent from kaons was 
measured in Ref. PQ, involving kaons produced at the primary vertex and kaons from 
bottom and charm decays. 



V Conclusion 

We have proposed a kinematic test to determine if the anomalous like-sign dimuon 
charge asymmetry recently reported by the DO Collaboration [TJ[2] is due to the decays 
of short-lived particles such as B mesons. Although charmed particles are also short- 
lived, they are not expected to contribute much to the like-sign signal. The test involves 
a reduction of the maximum impact parameter 60 to a value well above the vertex 
resolution but small compared with that characteristic of B decays. Although we have 
proposed bo = 100 yum, a dedicated simulation would probably be preferable in order to 
optimize this parameter. 
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